Analysis of microwave propagation in arrays of dielectric cylinders 
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We rigorously compute the propagation and scattering of microwaves by regular and defected ar- 
rays of dielectric cylinders in a uniform medium. Comparison with the previous experimental results 
is made, yielding agreements in the transmission coefficients for a certain range of frequencies. The 
results also show that localized states are possible for defected arrays. However, the experimentally 
claimed localized states are not observed. Moreover, comparison is made with previous theoretical 
results. The agreements and differences are pointed out. 
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Waves surround us. Direct human communication is 
mainly conveyed by acoustic waves, and is enriched by 
gestures which are passed into our eyes through optical 
waves. Nowadays electronic waves are also everywhere 
in our daily experiences such as audio & video systems, 
computers, and Nintendo games. Among many inter- 
esting properties pertinent to waves, the phenomenon of 
wave localization is probably most intriguing and remains 
as one of the unsolved mysteries of last century. 

The concept of wave localization was first proposed 
by Anderson nearly half a century ago, initially for elec- 
tronic systems Q, then extended to classical waves in- 
cluding microwaves to be addressed in this paper. It 
refers to situations that purely due to multiple scatter- 
ing of waves by disorders, transmitted waves are confined 
or frozen in space and the wave envelopes decay nearly 
exponentially along any direct ion 0. Since its inception, 
wave localization has attracted substantial efforts from 
all areas of expertise, signified by the great body of lit- 
eratures ranging from small scales such as electronic sys- 
tems, nano-structured materials, to large scales such as 
seismic and ocean waves, as reviewed and summarized 
in, for instance, Refs. 0,0,0,0,0,13,0,0, E3, EH 
In spite of the efforts, however, some fundamental issues 
concerning wave localization are still open |l3| . 

A few prominent problems in the investigation lie in 
the localization in two dimensional (2D) random media. 
First, although it has been conjectured, thought to be 
genuinely valid, by the scaling analysis 14| that all waves 
are localized in 2D random media for any given amount 
of disorders or impurities, the experimental observation 
of 2D localization is scarce. Even for the limited ex- 
perimental results that are available 15], the disparity 
OahJ still remains in the interpretation of the results. 
A main concern is whether the localization in 2D has in- 
deed been observed or not [16|. Some authors further 
pointed out that there still lacks a definite observation of 
2D wave localization, particularly for classical waves . 
In addition, how waves are localized in 2D and whether 
the localization features can be explained or fitted by the 
existing theory |19| are among a few important issues yet 
to be addressed. 

It would be too ambitious to answer all the unsolved 
issues dwelling on 2D localization in a single endeavor. 



However, a further look at some of previous experimen- 
tal observations of 2D localization, in the hope of clarify- 
ing some confusions, seems more realistic and imperative, 
and will therefore be the main theme of the present pa- 
per. 

Here we would like to concentrate on the experimen- 
tal observation of microwave propagation in 2D dielectric 
lattices, reported by McCall et al. [l5j. The reason is not 
only because the experiment is the first ever reported, but 
it has stimulated so many later investigations, and its im- 
portance could be inferred from the constantly increasing 
citation number. In |15| . the authors measured transmis- 
sion of microwaves through regular and defected arrays 
of dielectric cylinders. The transmission bands and pho- 
tonic band gaps for the regular arrays are found to be 
in excellent agreement with theoretical predictions. For 
the defected arrays, the authors suggested the observa- 
tion of a localized state. This observation, however, has 
been contrarily discussed and supported by two indepen- 
dent theoretical inspections 0,^3- ^ n this paper, we 
will use the standard multiple scattering approach first 
formulated by Twersky 20] to further examine the ex- 
perimental results. We will show that localized states 
are indeed possible for the defected arrays. However, the 
experimentally claimed localized states are not observed. 
Comparison will also be made with previous theoreti- 
cal results. Agreements and differences are found and 
discussed. Some contraries in the previous analysis are 
clarified. 

The systems considered here are from Ref. ^| and 
they are composed of almost ideal arrays of low-loss high 
dielectric constant cylinders. The cylinders are placed in 
parallel in an uniform medium to form a square lattice 
of lattice constant denoted as d. The ratio of dielectric 
constants between the medium and cylinders is 9. Two 
types of measurement are carried out in 15]. One is to 
measure the E-polarized microwave transmission across 
the samples which are shaped as a rectangular slab with 
thickness of 9d and width of 18d and the lattice con- 
stant is 1.27 cm. The other is to measure the intensity 
distribution at 11.2 GHz for a sample size of 7.6 cm x 
7.6 cm with a cylinder being removed from the center of 
the sample, i. e. defected array. However, the lattice ar- 
rangement for the second measurement differs from the 
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first in that the lattice constant is 1.59 cm; such an incon- 
sistency in the lattice arrangements has been discussed 
in 0. The authors in suggested that the localized 
states are observed for the defected array. 

In our computation, we will take the physical param- 
eters from |l5j. However, we will consider larger lattice 
sizes when appropriate to ensure the stability of the re- 
sults. Both lattice arrangements, i. e. d = 1.27 and 1.59 
cm corresponding to filling factors 0.4487 and 0.2863 re- 
spectively, will be considered. The cylinder radius is 0.48 
cm. In addition to computing the transmission across 
rectangular slabs, for comparison purpose we also com- 
pute the transmitted intensity when the source in placed 
inside a sample that takes a circular shape. The spa- 
tial distribution of energy will also be calculated to com- 
pare with the experiment. Moreover, the phase diagram 
method [2l| will be used to investigate a coherence phe- 
nomenon associated with localization. To be consistent 
with the experiment, we only consider the E-polarized 
microwaves. 



(a)d=1.27cm (b) d=1 .59 cm 




Band Structure Band Structure 

FIG. 1: The band structures of two 2D square lattices: (a) 
for d = 1.27 cm and (b) for d = 1.59 cm respectively. The 
dielectric constant is 9. The complete band gaps are apparent 
for both lattices. 



Figure ^ shows the band structures for the two 2D 
square lattices aforementioned. The results on the left 
panel reproduce that in (ljj . confirming our numerical 
codes. For later discussion, the band structures of lattice 
constant 1.59 cm are also plotted. Here we see that the 
two results are similar for frequencies below the second 
gap. There are three prominent complete bandgaps for d 
= 1.27, while there are only two complete bandgaps for 
d = 1.59. 

The transmission intensity is plotted as a function of 
frequencies in Fig. [3 The situations in (al), (a2) and 
(a3) are considered in 0. From this figure, we observe 
the following. First for d = 1.27 cm: (1) The results for 
situations considered by the experiment agree both qual- 
itatively and quantitatively well with the measurement 
for frequencies up to the beginning of the third bandgaps 
at about 14 GHz; (2) While we did not plot the case of 



propagation along the [21] direction, we plot the case of 
the transmission from a source located inside the sample, 
with a cylinder being deleted at the center (a4). The re- 
sult is in accordance with the case that the transmission 
is measured across the sample. In both cases, there is a 
transmission peak inside the valley at 11.27 GHz, slightly 
above the measured value, which is outside the complete 
gap, but within the gap along the [10] direction. The 
reason for the difference in the location and magnitude 
of the peak from the experiment is due to the resolution 
of plotting. (3) There are some significant differences be- 
tween the theory and experiment for frequencies starting 
from the third gap. For these frequencies, we observe 
clearly a well defined reduction in the transmission, well 
matching the band structures in Fig.^ an d the reduction 
is even deeper than the second gap. The experimental 
results, however, show some smeared reduction. Further 
calculations point to a few factors contributing to the 
disagreement, such as disorders, inhomogeneities in the 
radii of the cylinders, and the size of samples. 



(a)d=1.27cm (b)d=1.59cm 
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FIG. 2: The transmitted intensity is plotted as a function of 
frequency (in GHz) for (a) d = 1.27 cm shown in the left panel 
and (b) for d = 1.59 cm shown in the right panel respectively, 
(al & bl): Regular system with size of thickness 9d and width 
18d; the wave propagates in the [10] direction; both the source 
and the receiver are kept at one lattice constant away from 
the both ends of the sample. (a2 & b2): The same as (al & 
bl) except that the wave propagates in the [11] direction. (a3 
& b3): System with a single defect, i.e., the central cylinder is 
removed; the propagation is along the [10] direction; both the 
source and the receiver are kept at one lattice constant away 
from the both ends of the sample respectively. (a4 & b4): 
The source is moved into the center of the sample with the 
central cylinder being removed; the sample takes a circular 
shape of radius 9d. 



The results of the transmission through a slab with d = 
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1.59 cm are depicted on the right panel of Fig. [5] These 
diagrams show that the transmission is significantly pro- 
hibited in the second gap, consistent with the band struc- 
ture calculation. There is also a peak in the reduction 
regime. The frequency at which the energy distribution 
is measured is outside this forbidden regime. 

An important feature in Figs. (a3) and (a4) is the peak 
inside the second gap along the [10] direction but outside 
the complete bandgap, i. e. within the pseudo gap. Mc- 
Call et al. 0] attempted to explain this peak as the 
result of a localized state. To support this view, McCall 
et al. further measured the energy spatial-distribution 
at this peak and attempted, believed to be successful, to 
show that there is indeed a localized wave at this fre- 
quency. Unfortunately, this time they used the different 
lattice with d = 1.59 cm. The inconsistency has been 
pointed out in 16]. Meade et al. argued that the 
peak is due to the resonance of extended states, and the 
apparent localized wave shown by Fig. 4 in [l5j | is super- 
ficial and is due to the resonance. 

To resolve the disagreement between the measurement 
and the previous theoretical analysis, we show in Fig. |3| 
the 2D spatial distribution of energy at 11.27 GHz for 
the lattice of d = 1.27 cm and 1.59 cm with one cylinder 
being taken out at the center of the sample. We also plot 
the associated phase diagram. The source transmits at 
the center. 



(a) d=1.27cm 
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FIG. 3: Single defect: Phase and spatial energy distribution 
for frequency = 11.27 GHz with (a) d = 1.27 cm and (b) d 
= 1.59 cm respectively. Note that at this frequency there is 
a peak in the transmission for the system with d = 1.27 cm 
(c. f. Fig.|Ua3, a4)), while there is no such peak for d = 1.59 
cm (c. f. Fig. (b3, b4)). Hereafter, the trivial geometrical 
spreading effect is removed by the normalization with regard 
to the energy directly from the source. The black dots denote 
the cylinder. 



According to Refs. 0,H2, there will appear a coher- 
ence for the phase of the localized wave. We brief this 
feature here. The energy flow of EM waves is J ~ E x H. 
By invoking the Maxwell equations to relate the electrical 
and magnetic fields, we can derive that the time averaged 
energy flow is < J > t = — dtJ ~ \E\ 2 \79, where the 



electrical field is written as E — eE\E\e l9 , with eg denot- 
ing the direction, \E\ and 9 being the amplitude and the 
phase respectively. It is clear that when 9 is constant, at 
least by spatial domains, while \E\ ^ 0, the flow would 
come to a stop and the energy will be localized or stored 
in the space. The phase can be represented by a vector 
defined as v(r) = e x cos9 + e y sin 9, and the phase vectors 
at various space points can be drawn on the 2D coordi- 
nates. The coherence in the phases is a unique indication 
of localization. 

Fig. EUa) shows that the wave is indeed localized at 
11.27 GHz, as suggested in 15] for the case of d = 1.27 
cm. Here is shown that the energy is mostly confined 
within the lattice and the prescribed phase coherence pre- 
vails. The localization due to scattering by the cylinders 
is so strong that the magnitude of the intensity surpasses 
the energy directly from the source. However, this local- 
ized wave is not what was observed in the experiment. 
With the experimental setup of d = 1.59 cm, the ex- 
perimentally claimed localization is actually absent, as 
suggested by Fig. (b) where the phase vectors point to 
different directions indicating the flow of energy. Never- 
theless, the ripples in the energy distribution mimics that 
observed in the experiment. We note that for all cases, 
the sample size considered here is larger than that in the 
experiment. 




FIG. 4: The completely random array: The phase and energy 
diagram for frequency = 11.27 GHz. The black dots refer to 
the cylinders. 



To further confirm that localization is indeed possible 
at 11.27GHz, we also considered the case of a complete 
random array which is otherwise regular with d = 1.27 
cm. The results are shown in Fig. 0] Here the wave 
localization indeed appears. The energy is localized at 
the center of sample with phase vectors aligning along 
one direction. Due to the randomness, the peak value of 
the localized energy is reduced significantly. 

We further consider two other frequencies for the case 
of regular array d = 1.27 cm with one cylinder being re- 
moved: one is at 11.45 GHz, slightly above 11.27 GHz 
but still within the pseudo gap along the [10] direction, 
and the other one lies at 13 GHz between the second 
and third gaps. The results are illustrated in Fig. [5] 
They show that the localization is apparent for 11.45 
GHz while absent for 13 GHz. It is interesting to notice 
that the pointing direction of the phase vectors is mostly 
reversed from the case with the lower frequency, com- 
paring to Fig. 13 (a). It is also worth to noting that the 
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localization shown by both Figs. [3] and [5] occurs outside 
the complete band gap. These results will eventually give 
insight to the relation between gap and localization. The 
absence of localization in defected arrays may be due to 
the finite size of samples, in the view of previous stud- 
ies. But it might also hint that waves are not necessarily 
always localized in 2D, as suggested before |2l| . 



(a) f =1 1 .45 GHz 




(b) f=1 3.00 GHz 
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FIG. 5: Single defect for d=1.27 cm: Phase and energy dia- 
gram for (a) frequency = 11.45 GHz which is slightly above 
11.27 GHz at which we observed a peak in the transmission 
in Fig. [5] and for (b) frequency = 13 GHz which lies between 
the second and third gaps (see Fig. |Ja)). 
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FIG. 6: Transmitted intensity as a function of frequency for 
a disordered array of size 9 x 18. The waves propagates along 
the [10] direction. The randomness is twelve percent of ran- 
dom displacement of the cylinders within a circle around their 
original regular positions with regard to the lattice constant. 



Finally, to examine the difference between the theoret- 
ical and experimental transmission results when compar- 
ing the left panel of Fig. [21 and Fig. 3 in [l5|. we inves- 
tigate the effect of disorders. For brevity, we only show 
the effect of the positional disorder which is introduced 
by moving the cylinders around their regular lattice posi- 
tion to a certain extent. Here the disorder is measured as 
12% of random displacement of the cylinders within a cir- 
cle around their original regular positions with regard to 



the lattice constant. The result is shown in Fig. E] The 
result indeed shows that the disorder tends to degrade 
the transmission reduction in the third gap to a certain 
degree, while the disorder effect is not so significant for 
the second gap. 

In summary, in this paper we have shown a numeri- 
cal analysis of the previous experimental results on mi- 
crowave propagation through arrays of dielectric cylin- 
ders embedded in a uniform medium [l6| by the stan- 
dard multiple scattering method. Comparison has been 
also made with other theoretical results. The results re- 
veal that localized states are possible for defected arrays. 
However, the experimentally claimed localized states are 
not observed. The disorder effects are also presented. 
The present work is hopped to stimulate further experi- 
mental observations. 

Thanks to NSC and NCU for supports. 
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